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Abstract 

T— I , 

' Time-frequency transforms represent a signal as a mixture of its 

^ , time domain representation and its frequency domain representation. 

^ ' We present efficient algorithms for the quantum Zak transform and 

I quantum Weyl-Heisenberg transform. 

m ■ 

I ■ 

^ ! 1 Introduction 
O . 

f-*) . The Fourier transform is an operator that expresses a time-dependent sig- 

I nal as a sum (or integral) of periodic signals. In other words the Fourier 

^ i' transform changes a function of time s{t) into a function of frequency S{uj). 

' If a signal is a function of time it said to be in the "time domain" and 

^ ! if it is a function of frequency it is said to be in the "frequency domain". 

^ I For signals whose spectrum is changing in time, i.e. nonstationary signals, 

O^' sometimes the best description is a mixture of the time and frequency com- 

^ , ponents. Signal representations which mix the time and frequency domains 
are called, naturally enough, "time-frequency representations" and are often 



\^ • used to describe time- varying signals for which the pure frequency or Fourier 



representation is inadequate A familiar example of a time- frequency 

representation is a musical score, which describes when (time) certain notes 
(frequency) are to be played. 

Formally speaking for our present purposes, a quantum signal is sim- 
ply a quantum state \ip) where the Hilbert space is the group algebra C[G] 
of a finite abelian group G. The Quantum Fourier Transform (QFT) is 
central to the important quantum algorithms for factoring and discrete log- 
arithm. Mathematically speaking, the Quantum Fourier Transform is a 
linear operator on the Hilbert Space C[G] which is a change of basis from 
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the basis of group elements {\gi), I^igi)} to the basis of characters of G, 

{|Xl),IX2>,.--,|X|G|>}- 

We present efficient algorithms for quantum versions of the Zak and 
Weyl-Heisenberg transforms. Both these time-frequency tranforms can be 
seen as generalizations of Fourier transforms and the quantum algorithms 
make heavy use of the Quantum Fourier Transform. We follow the theory 
and notation of Q and recommend this book as background to this material. 

2 Zak Transforms 
2.1 Background 

Let A be a finite, abelian group, A* the group of characters of A (note: 
in this paper * does not mean conjugation), B < A a subgroup of A, 
= {a* G A* : a*{b) = l,b £ B] the dual to B, and / e C[A], the 
group algebra of A. Define 

Z{B)f G C[A X A*] 

by the formula 

Z{B)f{a,a*) = Y,f{a + b)^. 

beB 

F = Z{B)f is called the Zak transform of / over B. A simple calculation 
shows that F{a + b,a* + b^) = a*{b)F{a,a*) where b £ B and 6* € 5*. 
Therefore F is determined by its values on a set of coset representatives of 
B ^ B^ m. A X A* and thus conceptually we may think of F as a function 
on T where T is a set of coset representatives. Since 

I B X Bn, I I i? 1 1 B^ I 

we have the same number of degrees of freedom with which we started. 
Notice that if B contains only the identity, i.e. is the trivial subgroup, then 
Z{B)f{a,a*) = f{a) and is basically the identity map. Also notice that 
if B = A then Z{A)f{0,a*) = {a*\f) and therefore Z{A)f is basically the 
Fourier transform of /. So the Zak transform mediates between the time 
domain and frequency domain depending on the subgroup B. 

Consider the function /(ao) = 5{x — gq) which is 1 on ao and 
otherwise. Applying the above formula for the Zak transform yields 
F{a,a*) = a*{a — ao) for a G aQ + B,a* £ A* and otherwise. But since 
F is determined by its values on a set of coset representatives of B x B^ in 
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Ax A* let us introduce such a set of representatives T = T1XT2 = {{xi,aj)} 
where Ti = {xi} is a set of coset representatives of -B in ^ and T2 = {a^} is 
a set of coset representatives of B^, in A* . Bearing in mind the above trans- 
formation of a delta function, we now offer our definition of the Quantum 
Zak Transform (QZT) (with respect to T) by 



where Xa G Ti is the coset representative of a. Now notice that Xa — a £ B. 
Therefore a| is restricted to B and therefore can be considered to be a 
character of B, i.e. an element of B*, and this restriction is independent of 
the choice of coset representative, i.e. it is natural or canonical. Therefore 
an equivalent formulation of the QZT is given by 



a)^— ^ b*{xa-a)\xa)\b*). 



The only difference in these two formulations is in the interpretation of the 
observed content of the second register. 



2.2 The Quantum Algorithm 

We now show that the QZT is efficiently implementable. Define P{B) to be 
the transform 

P{B)\a) = \xa)\xa-a) 

which decomposes a into its coset representative and the corresponding ele- 
ment of -B. P is clearly unitary and efficiently implementable. After applying 
P(B) we apply the Quantum Fourier Transform (over the group B, denoted 
Fb) to the second register. This results in the state 

^ b%Xa-a)\Xa)r). 

Therefore the QZT is simply Z{B) = {10 Fb) o P{B). 



3 Weyl-Heisenberg Transforms 
3.1 Background 

Define g(^x,x*){o) = ff(^ ~ x)x*{a) to be the time-frequency translate of g by 
(x, X*) where g G C[A\. We will use time-frequency translates to form 
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orthonormal bases so we also require \g\ = 1. Let A = B x and 
((7, A) = {g(^x,x*) '■ ib,b^) £ A}. We call (5, A) a W-H system over A with 
window g. A basic result (Q, Theorem 12.1 corrected version) is that ((7, A) 
is an orthonormal basis of C[A\ if and only if for all (a,o*) S A x A* we 

have |G(a,a*)| = -^/j^ where the Zak tranform is taken over B. Because 
von Neumann measurements must be unitary we will restrict our attention 
to window functions g which satisfy this constraint. Utilizing POVMs one 
could consider implementing nonorthonormal W-H systems but we will not 
address this in this note. This orthogonality constraint together with the 
earlier observation that G is determined by its values on a set of coset rep- 
resentatives of B X B^ in Ax A* implies that orthonormal W-H systems are 
in bijective correspondence with the set of all |j4|-tuples of complex num- 



bers with modulus y]^- Iii this note we will restrict the W-H systems 
under consideration by assuming that for each (a,a*) G A x A* the phase 
of G{a,a*) is a rational fraction of 2tt which we can compute in polynomial 
time. Whether or not this last assumption is excessively restrictive would 
depend on the intended application. Notice that if g is the constant function 
g = J and A = {0} x A* then ((7, A) is the (normalized) Fourier basis, 

G(a,a*) = ^1^^ and this restriction holds trivially. 

We define the Quantum Weyl-Heisenberg Transform (QWHT) by 

(fe,b.)GA 

In other words, the QWHT expresses 1-0) in the orthonormal basis of time- 
frequency translates of the window function. 

3.2 The Quantum Algorithm 

Let 

(b,6*)eA 

i.e. the a's are the coeffients of the WH-expansion of /. Define 

P{a,a*)= aib, b^)b^{a)a*{b). 

Notice that P is A-periodic and that the a's are, by definition, the Fourier 
coefficents (over A x A*) of P. A fundamental result (0], Theorem 7.5) 
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states that F = GP. This result suggests an algorithm for computing the 
WH-coefficients of /, namely compute the Fourier coefficients of P = ^. 

Define ^{g) to be the unitary transformation which acts on the Hilbert 
space C[T] (recall T is the set of coset representatives B x B^, in A x A*) 
by 

Since the phase of G{xi,a*) is, by assumption, a rational fraction of 27r 
computable in polynomial time we may efficiently implement ^{g) by the 
phase kickback technique described in Finally in order to complete our 
description of the algorithm, we must assume that we are given an explicit 
isomorphism between A and A*. These groups are isomorphic, though not 
canonically so. Therefore in any computational situation we provide an 
explicit isomorphism by choosing an explicit computational representation 
of the groups A and A*. This isomorphism induces explicit isomorphisms 
between B and B* and between the factor group A/B and B^. We will see 
shortly how we will employ these three interrelated isomorphisms. We will 
highlight this interrelation, and abuse notation, by using the symbol (j) to 
refer to all three of these isomorphisms, allowing for context to make the 
usage clear. As in the case of the Zak transformtion, these isomorphisms 
are simply reinterpretations of the contents of the registers. 

Our QWHT is the sequence Fb,xB ° ^{g) ° ^{B)- Let us see how this 
unitary transformation acts on \a). We have 

Z{B)\a) = — ^ ^ a*{xa - a)\xa)\a*) 



and then after applying ^{g) we obtain: 

which by the fundamental result discussed above equals: 



y\B\ I,* 

where we are now considering the contents of the second register to be an 
element of B* . We now utilize our explicit isomorphisms to reinterpret the 
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contents of the first register as an element of and the contents of the 
second register as an element of B: 



B 



By applying tlic final transformation in the sequence Fb^xB we obtain our 
desired expansion: 

E i'^\9(bM)\b)\b*)- 
(6,6*)eA 
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